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Abstract 

In a roughening process, the growth exponent (3 describes how the rough- 
ness w grows with the time t: w ^ t^. We determine the exponent /3 of a 
growth process generated by the spatiotemporal patterns of the one dimen- 
sional Domany-Kinzel cellular automaton. The values obtained for (3 shows 
a cusp at the frozen/active transition which permits determination of the 
transition line. The /? value at the transition depends on the scheme used: 
symmetric (/3 ~ 0.83) or non-symmetric (/3 ~ 0.61). Using damage spreading 
ideas, we also determine the active/chaotic transition line; this line depends 
on how the replicas are updated. 
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I. INTRODUCTION 



The one dimensional Domany-Kinzel cellular automaton (DKCA) is a totally discrete 
systems - temporally, spatially and in numbers of states - with several applications in physics, 
chemistry, biology, computer science, etc. The DKCA phase diagram was originally proposed 
by Domany and Kinzel |]T| who showed the existence of two phases: an active phase and 
a frozen one. A more detailed study, using numerical simulation, was present by Martins 
et al 0], where a new phase in the active region - a chaotic phase - was discovered using 
the damage spreading technique. Further, Zebende and Penna used the gradient method 
to determine the phase boundaries with high precision. Recently, Tome explored some 
details of the joint evolution of two DKCA. She considered the problem of simultaneous 
updating of two replicas using pseudo-random numbers; two prescriptions were presented 
for the joint evolution: one used by Martins et al 0] and another introduced by Kohring 
and Schereckenberg . Hinrichsen et al [^] discovered a third phase in the DKCA diagram, 
using a thorough analysis of the damage spreading technique to spht the active phase into 
three different regions: a chaotic region, where the damage spreads for every member of this 
family of dynamic procedures; an active region, where the damage heals for every member 
of this family, and another active region, where the damage spreads for a subset of the 
possible dynamic procedures and heals for the others. This new phase was obtained with a 
prescription that updates the replicas using the minimal correlations. 

In 1997, de Sales et al showed that the roughness exponent a can be used to classify 
the deterministic cellular automata (CA) described by Wolfram [§. More recently [^, they 
showed that this exponent also can be used to detect the frozen/active transition in DKCA 
directly from the automata, without reference to order parameters or response functions. 
This method also can be used to detect phase transitions in another kinds of models, such as 
the Potts model Beyond the roughness exponent a, the growth exponent (3 is another 
critical exponent used to describe various roughening processes in the surface growth context 

ini,n:2i. 
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In this work, we introduce the growth exponent method to identify phase transitions. 
We apply this method to the one dimensional DKCA to detect the phase transitions directly 
from the automata and build the DKCA phase diagram. To obtain self-affine rough profiles, 
we use the accumulation method to perform a mapping to the solid-on-solid (SOS) model, 
and measure the time evolution of the roughness of the interface to obtain the exponent f3. 
In the frozen/active transition, the exponent (3 has a maximum, and two schemes are used 
to update the system [Q: a symmetric scheme, that corresponds to a triangular lattice, 
and a nonsymmetric scheme. In the symmetric scheme we obtain f] ~ 0.83, consistent with 
the directed percolation (DP) prediction; but in the nonsymmetric scheme, (3 ~ 0.61, an 
unusual value that was not expected. It is expected |T^, rather, to find the same value for 
/3 in the two schemes. 

Very recently, a quite similar method was used by Lauritsen and Alava |T^, to study 



the Edwards- Wilkinson equation with columnar noise, and by Dickman |T^, to study the 
contact process. 

This studies do not consider damage spreading, and cannot detect the chaotic/non- 
chaotic transition. To evidence this transition, we use the damage spreading ideas []TB 



the difference between two replicas, which evolve with same dynamics, is used to perform 
the same SOS mapping as mentioned above, via the accumulation method. We use three 
different prescriptions for the simultaneous updating of the replicas and obtain three different 
transition lines. 

In Section II, we introduce the growth exponent method and show the results obtained 
for the frozen/active transition. The damage spreading ideas and the chaotic/non-chaotic 
transition are presented in the Section III. Finally, we let to the Section IV our conclusions 
and acknowledgements. 
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II. FROZEN-ACTIVE TRANSITION 



The DKCA consists of a linear chain of L sites {i = 1, 2, L), with periodic boundary 
conditions, where each site has two possible states a = 0, 1 (frozen, active). The state of the 
system at time t is given by the set {cxj}. In contrast to the deterministic CA studied by 
Wolfram the DKCA is probabilistic, in the sense that the rules that update the system 
are given by conditional probabilities, P((Ji_i(t), (7j+i(t)|(Tj(t+l)) (in the symmetric scheme). 
That is, the state of given site in time {t + 1) depends in a probabilistic fashion upon the 
state of the two nearest neighbours at time t. In the DKCA, the conditional probabilities 
in the isotropic case are 

P(0,1|1) = P(1,0|1) =pi, P(l,l|l) =p2, P(0,0|1) = 0. 

Obviously, 

P((Ti„l,(Ti+i|0) = 1 - P((Ti_i,(Ti+i|l) . 

In the nonsymmetric scheme, the conditional probabilities are given by P(crj_i(t), ai(t)\ai{t+ 
1)), and P((jj_i, cTjlO) = 1 — P(o"i_i, (Tj|l) . The values of the conditional probabilities are 
the same in both schemes. In Figure 1 we reproduce a representation of these schemes. 

Depending on the values of the parameters {pi,P2), the asymptotic {t — > oo) state of the 
system is either a frozen state, with all sites in state 0, or has a finite fraction of sites with 
value 1, the active state. This is a second order phase transition, characterized by universal 
critical exponents 

To study the phase diagram of the DKCA, we use the accumulation method, introduced 
by Sales et al. 0, to obtain profiles in analogy with solid-on-solid (SOS) models in (1 + 1) 
dimensions. This method consists in accumulating (or summing) all the values assumed by 
the variables aiit) during a given number t of successive time steps 

h,{t) = iKr) . (1) 

r=0 
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The differences between the schemes become exphcit at this point. In the Figure 2, 
we can observe several profiles at criticality in both schemes, generated after we apply 
the accumulation method. It is evident that different schemes lead to completely different 
profiles. 

Thus, we obtain growth processes, the nature of whose correlations can be investigated 
through the analysis of the roughness w{L,t) The roughness is defined by 

L 



w'iL,t) = ^j:{Ht)-Kt)y , (2) 

i=l 

where h{t) is the mean value of hi{t) at time t. In fact, in order to consider the initial 



roughness of the profiles, we work with the fluctuation in roughness [17 



Sw{L, t) = ^w^{L,t) -w'^{L, 0) . (3) 

We expect that the behavior of Sw{L,t) has the form 

Swmt)^L-f^-^^ , (4) 

where f{u) is a universal scaling function, a is the roughness exponent, z = is the 

dynamic exponent and (3 is the growth exponent. The function f{u) = constant at large 
times (t ^ L^) and f{u) ~ at short times {t <^ L^). So, at short times, we expect that 
6w{L, t) ~ t^ and can measure /5 calculating the slope of the log — log plot of 6w{L, t) versus 
t. The growth exponent denotes how the profile roughness grows with time: (3 = 1/2 means 
that the profile is not correlated, and is analogous to that generated by random deposition 



rT|]; if /9 > 1/2, the profile tends to grow more at the tips, which causes the roughness 
increases faster, in contrast to f3 < 1/2, where the valleys grows quickly and make the 
roughness increase more slowly. 

Typical results for the evolution of the roughness are showed in Figure 3. These results 
correspond to averages over 50 random initial configurations taken after 100000 time steps 
in DKCA containing L = 10000 sites. Each curve in this graph (101 points) takes approx- 
imately one hour of CPU time on a Digital SOOaw workstation. We can observe that the 



roughness reaches the steady state in the frozen phase and grows indefinitely in the active 
phase. The values of (3 at the transition are showed in Figure 4, for both schemes. The ex- 
ponent (3 is measured over more than three decades (10 <t< 100000 and 0.1 < 5w < 1000). 
Note that (3 shows a maximum at the transition and tends to the value (3 = 1/2 quickly 
after the transition, remaining at this value until pi = 1. The value of the (3 exponent at the 
transition depends of the scheme utilized in the DKCA. Nagy et al describe two schemes: a 
symmetric one, where the even (odd) sites are updated at even (odd) times; and a nonsym- 
metric one, where all sites are updated at each time step, but the neighbours of site («, t + 1) 
are [i — l,t) and (i,t). The symmetric scheme has been used in most previos studies of 
the model. In the symmetric scheme, we find j3 = 0.81(2), compatible with the universality 
class of directed percolation (DP) {f3 ~ 0.84) [|15|. In the nonsymmetric scheme, we find 
/3 = 0.61(2), clealy different from the DP value. 

A finite size scaling analysis was made for the exponent (3 and shows that the width of 
the peak vanishes when the size of system goes to infinity. The value of the exponent (3 
when the system size goes to infinity approaches (3 ~ 0.83(2) at the frozen/active transition 
in the symmetric scheme, and is valid for all values of p2 7^ 1- 

On the line p2 = 1, the system is mapped in the two-dimensional Ising model []TB|, and 



the f3 value is significant greater in the two schemes: (3 = 0.99(1) in the symmetric, and 
(3 = 0.75(1) in the nonsymmetric. The value in the symmetric scheme agrees with the 
literature 0, which predicts that at all points on phase boundary, except the line p2 = 1, 
are characterized by directed percolation (DP) exponents. On the line p2 = 1, the model 
has been solved exactly, and can be mapped in the bi-dimensional Ising model universality 
class, which leads to the value /3 = 1 |[l|]. 

The value (3 > 1/2 denotes the trend of the system to grow faster at the tips which can be 
understood as a preservation of active sites. At the transition, few sites remain active, which 
causes hi{t) to grow only near active sites, contributing to an increase in the roughness. In 
the active phase, many sites are active, but uncorrelated, randomly increasing the height 
hi{t), so that (3 = 1/2. 
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The roughness exponent is defined only in the frozen phase because, in the active phase, 
the steady state is not reached. We calculate the Hurst exponent, if, very close to the 
transition and find, in the symmetric scheme, H = 0.61(2). As noted by de Sales et a/ , 
there is a maximum in the H exponent, marking the phase transition. In the nonsymmetric 
scheme, the value of H at the transition is much smaller, a ~ 0.25, and the maximum is not 
very clear. 

III. ACTIVE-CHAOTIC TRANSITION 

Martins et al [0] used the damage spreading technique to show that the active phase of 
DKCA can be split into two phases, chaotic and non-chaotic. The order parameter of this 
transition is the difference between two replicas with slightly different initial configurations. 
They let the system evolve until it attains equilibrium, and then a replica of the automaton 
is created with some sites altered (damage). So the two replicas, one with states ai{t) the 
another with states Qi{t), evolve with the same dynamics, and the difference between the 
automata 

is measured. The fraction of sites in replica system that differ from their counterpart in the 
original system is called the Hamming distance, defined as 

DH{t) = WUt) ■ 

^ i 

The stationary Hamming distance is null in the non-chaotic phase and positive in the chaotic 
phase. 

To obtain the chaotic/non-chaotic boundary, we use a slightly different method, where 
the difference between the two automata is used to perform the same mapping to a SOS 
model as we did in the accumulation method 

K{t) = jzT,{T) . (5) 

r=0 
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Thus, the profile generated by the difference of the two replicas behaves exactly as the 
profiles generated in the frozen/active boundary: the roughness reaches a stationary value 
in the non-chaotic phase and grows indefinitely in the chaotic phase. This behavior can be 
understood if we note that the difference between the replicas vanishes in the non-chaotic 
phase, which implies no contribution to the height hi{t), and is positive in chaotic phase, 
which implies in a persistent contribution to the height. The (3 exponent again passes 
through a maximum at the chaotic/non-chaotic transition, and its value depends on the 
scheme utilized. 

Figure 5 shows (3 for the chaotic/non-chaotic transitions, in the symmetric scheme. To 
obtain this figure, we use an automaton with L = 10000 and let it evolve for 10^ time 
steps, at which time we create a replica of the system with an "initial damage" , by flipping 
a fraction of the sites (~ 10%). Then the replicas evolve with the same dynamics during 
the 10^ time steps, and the difference between them is measured as a function of time. An 
average over 50 realizations of initial damage was used. 

We can locate the chaotic/non-chaotic phase transition, with the damage method, either 
waiting or not waiting for the original system to reach the steady state (first 10^ time 
steps). The only difference we note is that the value for the exponent (3, at the frozen/active 
transition, is more precise and the cusp is more pronounced if we wait the original to reach 
the steady state. For the chaotic/non-chaotic transition we obtain the same values with the 
two procedures, considering the statistical fluctuations. 

At this point we have to emphasize the question of the dynamics of joint evolution of two 
OA's. Tome studied this joint evolution and showed that the CA's can evolve following 
two different prescriptions: prescription A, used by Martins et al. p[, which corresponds 
to updating both automata always using the same random number, and prescription B, 
introduced by Kohring and Schereckenberg [Q, which implies that one must sometimes use 
two different random numbers [zi and Z2) to update the original and the replica. The later 
occurs when we have (cTj-i + cTi+i) = 1 and + Qi+i) = 2, or vice-versa. We make 
simulations using these two prescriptions and verify that there are significant differences in 
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the chaotic/non-chaotic boundary of the DKCA p2 —pi phase diagram, as showed in Figure 



6. These differences were discovered by Bagnoh [^], who studied the damage spreading 
transition in the DKCA using a mean-field approximation; this work confirms his prediction 
numerically. The frozen/active transition line and the chaotic/non-chaotic transition line 
obtained with prescription A are very close to the the phase diagram found by Zebende and 
Penna ||. 

In this phase diagram we also present a third prescription, C, in which three random 
numbers, Zi, Zi and 2:3, are used to update the system. Defining 

IJi = ai-i + 2 * ai+i 

and 

Vi = Qi-i + 2* Qi+i , 

we have the following cases: 

• if Ui = Vi ^ we use the same random number (zi) to update the original and the replica; 

• if f/j = 1 and = 2 (or vice-versa) — ^ we use zi for the original and Z2 for the replica; 

• if f/j = 1 and = 3 (or vice- versa) we use zi for the original and z^ for the replica. 

• if f/j = 2 and Vi = 3 (or vice- versa) we use Z2 for the original and z^ for the replica. 
To this prescription, the boundary is slighting above that of prescription B. 

The difference between the phase boundaries appears to be due the different prescrip- 
tions for updating the systems in the damage spreading technique employed to detect the 
chaotic/non-chaotic transition: prescription A corresponds to maximal correlations between 
the random numbers; prescription B, to lower correlations, and prescription C to minimal 
correlations P,p!9|. We perform simulations using the third prescription in an attempt to 
reveal the third phase transition, reported by Hinrichsen et a/ 0, but are unable to detect 
that transition using prescription C, which is the prescription with minimal correlations. 
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IV. CONCLUSIONS 



In this work we propose a new method to obtain the phase diagram of the DKCA using 
the growth exponent /?, and detect the frozen/active boundary. At the transition, the values 
of the exponent (3 depend on the scheme used to update the system: (3 = 0.83(2) for the 
symmetric scheme, and (3 = 0.61(2) to the nonsymmetric scheme. Next, we extend this 
method, using damage spreading, to obtain the chaotic/non-chaotic boundary. Finally, we 
study three different prescriptions for the joint evolution of two DKCA's and construct the 
phase diagram, showed in Figure 6. 

The advantage of this method to determine the phase diagram of the DKCA is that we 
do not need to wait for the system to reach the steady state, as in the methods used by 
Martins et al and Zebende and Penna [0], thereby economizing computation time. In 
addition, the growth exponent method can detect the chaotic/non-chaotic boundary much 
more clearly than the usual Hamming distance, which presents large fluctuations at the 
transition. This method can also be employed to detect phase transitions in other models 
where the accumulation process can be used [jlO| . 
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Figure captions 

Figure 1 Schemes used to update the systems. 

Figure 2 Evolution of profiles generated by the accumulation method. Above: profiles 
generated in the symmetric scheme on a lattice with L — 1000. Below: profiles generated 
in the nonsymmetric scheme, L — 1000. Both figures have the same initial profile and same 
sequence of random noise, and the profiles were taken at the same instants of time in both 
schemes. Each profile was taken after 1000 time steps. 

Figure 3 Evolution of the fluctuation in roughness 5w{L,t) with time i in a log-log plot, 
for p2 = 0.95 and five different values of pi. We use L = 10000 and 50 samples. Note that 
the frozen/active transition occurs when the roughness grows indefinitely (pi = 0.595 - filled 
symbols) . 

Figure 4 Evolution of growth exponent /3 as a function of pi, in the frozen/active transition, 
in the two updating schemes; the system has L — 10000. Five different values of p2 are shown. 
The maximum of /3 indicates the transition. The symmetric scheme is represented by filled 
symbols and the nonsymmetric scheme by open symbols. The lines are to guide the eye. 
Figure 5 Evolution of growth exponent /3 as a function of pi, using damage spreading ideas 
to locate the chaotic/non-chaotic transition. The system has L — 10000. Four different 
values of p2 are showed, for prescription A. 

Figure 6 DKCA pi — p2 phase diagram obtained via the growth exponent method. Note that 
the three prescriptions yield different boundaries for the chaotic/non-chaotic transition. At 

P2 = we have pic = 0.809(1) and at pi = 1 we have p2c = 0.3135(10), with all boundaries 
meeting at these points. 
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